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Abstract 

The theory of quantum propagator and time-dependent integrals of motion 
in quantum optics is reviewed as well as the properties of Wigner function, Q- 
function, and coherent state representation. Propagators and wave functions 
of a free particle, harmonic oscillator, and the oscillator with varying frequency 
are studied using time-dependent linear in position and momentum integrals 
of motion. Such nonclassical states of light (of quantum systems) as squeezed 
states, correlated states, even and odd coherent states (Schrodinger cat states) 
are considered. Photon distribution functions of Schrodinger cat male and fe- 
male states are given, and the photon distribution function of squeezed vacuum 
is derived using the theory of the oscillator with varying parameters. Properties 
of multivariable Hermite polynomials used for the description of the multimode 
squeezed and correlated light and polymode Schrodinger cats are studied. 

1 Introduction 

In the lectures the recent results relating time-dependent integrals of motion 
and Green function of quantum system (Feynman path integral) are demonstrated. 
Obtaining of the quantum propagator of a quantum system with quadratic Hamilto- 
nians is reduced to solving the system of the eigenvalue equations for the linear in 
positions and momenta integrals of motion. In spite that in the Schrodinger equation 
for the propagator the time and space variables are not separated the existence of the 
time-dependent invariants gives the possibility to treat this problem as the problem 
with separable variables. 

The integral of motion which is quadratic in position and momentum was found 
for classical oscillator with time-dependent frequency a long time ago by Ermakov 
0. Two time-dependent integrals of motion which are linear forms in position and 
momentum for the classical and quantum oscillator with time-dependent frequency 
were found in 0]; for a charge moving in varying in time uniform magnetic field, this 
was done in |^ (in the absense of the electric field), [^j, 0] (nonstationary magnetic 
field with the "circular" gauge of the vector potential plus uniform nonstationary 



electric field, and P (for the Landau gauge of the time-dependent vector potential plus 
nonstationary electric field). For the multimode nonstationary oscillatory systems, 
such new integrals of motion, both of Ermakov's type (quadratic in positions and 
momenta) and linear in position and momenta, generalizing the results of @] were 
constructed in 0. Recently the discussed time-dependent invariants were obtained 
using Noether's theorem in |]10|-|]12 . 



In the lectures the coherent states of the oscillator are discussed and these states 
are used to construct Q-function of the oscillator at zero and finite temperatures. 
The nonclassical states are discussed in details. The Wigner function of multimode 
squeezed light is studied using such special functions as multivariable Hermite poly- 
nomials. The theory of parametric oscillator is appropriate to consider the problem 
of creation of photons from vacuum in a resonator with moving walls (with moving 
mirrors) which is the phenomenon based on the existence of Casimir forces (so-called 
nonstationary Casimir effect). The resonator with moving boundaries (moving mir- 
rors, media with time-dependent refractive index) produces also effect of squeezing in 
the light quadratures. In the high energy physics very fast particle collisions may pro- 
duce new types of states of boson fields (pions, for example) which are squeezed and 
correlated states studied in quantum optics but almost unknown in particle physics, 
both theoretically and experimentally. 

2 States of Quantum Systems 

In the Schrodinger representation, the state of one-dimensional quantum system is 
described by a complex wave function (in coordinate representation) 

iP{q,t) =<q\ij,t>, (1) 

where {ip, t > is the time-dependent state vector in Hilbert space of quantum system 
states, and |g > is the eigenstate of the position operator. Such states are called 
pure states. 

The mixed state of the quantum system is described by the density matrix (a 
function of two variables q,q', and time t) 

p{q,q',t)=<q\p{t)\q'>, (2) 

where p(t) is the time-dependent density operator satisfying the conditions of her- 
mit icity 

pm = pit), (3) 

normalization 

Tr Pit) = 1, (4) 

and nonnegativity 

Pit) > 0. (5) 



For pure state, the density operator which is projector on the state \ip,t >, i.e. 

p(t) = \^,t><ij,t\ (6) 

satisfies the extra condition 

p\t)=p{t), (7) 

and, consequently, 

lx = Trp\t) = l. (8) 

For the system with a Hamiltonian operator H{t), the state vector satisfies the 
Schrodinger equation ( h is Planck constant) 

tn^^\^,t>=H{t)\^,t>. (9) 

The density operator (||) obeys the evolution equation 

^ + ^[^W,pW] = o. (10) 

All the information on the quantum system is contained in the density operator whose 
matrix elements may be taken in any appropriate basis. 

3 Evolution Operator 

The evolution of the state vector {ip, t > is determined by an evolution operator U{t) 
which is defined by the equality 

\^,t>=U{t)\^,0> . (11) 
This operator obeys the Schrodinger equation 

= Hit)U{t), (12) 

with the obvious initial condition 

f/(0) = i, (13) 

where 1 is the unity operator. The matrix element of the evolution operator in coor- 
dinate representation is called Green function (or quantum propagator, or Feynman 
amplitude considered as path integral) 

G{q,q',t)=<q\Uit)\q'> . (14) 



The Green function satisfies the Schrodinger equation 



^^ dG{q^q',t) _ ^ thSmq - q'), (15) 

with the initial condition 

Giq,q',0) = 6iq-q'). (16) 

The propagator in momentum representation is connected with the Green function 
by the integral transform 

G{p,p',t) =<p\U{t)\p' >= J <p\q>< q\U{t)\q' >< q'\p > dq dq\ (17) 

which is the Fourier transform since 

1 ipq 

and we used the completeness condition of the position state vectors 

\q X q\ dq = 1. (19) 
Below we will use the propagator and density matrix in other representations, too. 

4 Integrals of Motion 

The integral of motion or invariant operator I[t) is defined by the equality 

j^<tlj,t\l{t)\tlj,t>=0. (20) 

The sufficient condition for the operator to be the invariant is the equality to zero of 
the full time derivative of this operator, 

Comparing Eq. (PT| ) with Eq. ( PUj) we see that the density operator p{t) is the 
quantum integral of motion for systems with Hamiltonian H{t). 

There exist several theorems on the properties of integrals of motion 

i) A function of integrals of motion F{I{t)) is the integral of motion. 

ii) Eigenvalues of integrals of motion do not depend on time. 

iii) Given a solution to Schrodinger equation |^/'(t) > and integral of motion I{t). 
Then the vector |0, t >= I(t)\ilj{t) > is solution to Schrodinger equation. 



These properties are easily proved if one uses necessary and sufficient conditions 
for the operator I{t) to be the integral of motion which are expressed by the relation 
of the quantum invariant and the evolution operator of the system U{t) of the form 

I{t) = U{t)I{0)U-\t). (22) 

Systems with hermitian Hamiltonians H^{t) = H{t) have unitary evolution operators 

U{t)^ = U-\t). (23) 

Observables I nit) in the Heisenberg representation are related to the evolution op- 
erator by the formula 

lH{t) = t/t(t)/(0)t/(t). (24) 



From Eqs. (|22D-(p^, the relation of integrals of motion I(t) to observables in the 
Heisenberg representation is [0 

I{t) = U\t)lH{t)U-\t). (25) 

For stationary systems with time-independent Hamiltonians, this relation is reduced 
to the equality 

lit) = Ini-t). (26) 

There are integrals of motion which at the initial time t = coincide with the position 
and momentum operators. These invariants are 

Poit) = Uit)pU-'it), (27) 

and 

qo{t) = Uit)qU-\t). (28) 



5 Invariants and Propagator 

From relations ([27|), ( pSj) and the rule of the multiplication of matrices, it follows the 
system of new equations for the quantum propagator |l| 

Po{t)gG{(l,Q',t) = ^f^-^G{q,q',t), (29) 

qo{t)M(l,<l'^t) = q'G{q,q',t). (30) 

In the left hand side of equalities subindices q mean that operators Po{t) and qo{t) 
act on the argument q of the Green function. The system of Eqs. (0), (|29|) , ( pOD 
is selfconsistent system, and the coordinate dependence of the quantum propagator 
may be found from Eqs. (p9|), (^). The time-dependent phase factor of the Green 
function is determined by Eq. ([T5|). Analogous systems of equations for the quantum 
propagator may be written in any representation. 



6 Free Particle 

The Hamiltonian for a free particle with mass m is time-independent 

H=^. (31) 

The momentum operator is known to be the integral of motion for the free motion. 
So we have 

Po{t)=p. (32) 
It is easy to check that the time-dependent operator 

qo{t) = q~-t (33) 
m 



is the integral of motion since its full time derivative equals to zero. Equality ( ^6]) 
means that the position and momentum operators of the free particle in the Heisen- 
berg representation are 

PH{t)=p, (34) 

qH{t)=q + -t. (35) 
m 

If one introduces the vector-operator Q with components Qa-, a = 1, 2, such that 
Qi = Pi Q2 = q, Eqs- (HI)) (HI) Hiay be written in the matrix form 

Qo(t) = A(t)Q, (36) 

where the 2x2-matrix A(t) has the form 

= ( -t/,„ 1 ) ■ P7) 
One can check that the commutator 

[Qa, Q/3] = So/3, (38) 

where the matrix S has the form 

E = -tn(^^^ J), (39) 

is not changed for integrals of motion 

[Q0a{t),Q0f3{t)]=^ap. (40) 

It means that the matrix A(t) is the real symplectic matrix, i.e. 

A(t)SA*"(t) = S (41) 



belonging to the Lie group Sp{2, R). 

It is interesting to note that the time-dependent operator quadratic in position 
and momentum 



Kit) = g2 + ^ _ l^pq + (42) 

is the integral of motion for the free particle. This operator is obtained using the 
relation 

m = qlit) (43) 
analogous to the relation of the energy to the invariant momentum po{t) ( p2D 

According to Eqs. (p9[), (|30D and (0), (|3^) the free particle propagator must satisfy 
the system of equations 

' dq dq' ' ^ ' 

ni ' +\ ^ 9G{q,q',t) , 

qG{q,q,t) + = qG{q,q,t). (45) 

m oq 



Equation (Q) means that 

G{q,q',t) = g{x,t), x = q - q' , (46) 
where g{x, t) is an arbitrary function. Equation (H3|) may be rewritten in the form 



ihtdqixA) , 

xg{x,t)^ ^^ = 0- 47 

m ox 



Then we obtain the solution for this equation in Gaussian form 



2 

^7(x,t) = iV(t)exp^. (48) 

The function N{t) is obtained from Eq. (p!5|) 

N{t) = ^. (49) 

The normalization constant G can be found from the initial condition for the prop- 
agator ([T6|) using the formula 

1 

lini{— ;=exp } = (5(x), (50) 



which gives the free particle propagator 



7 Harmonic Oscillator 

Let us study by methods of the integrals of motion the harmonic oscillator with the 
time-independent Hamiltonian 

i7 = ^ + ^. (52) 

One can check that there exist two invariants linear in position and momentum 

Po(t) = cosLut p + musmut q (53) 

and ^ 

Qoit) = smut p + cosut q. (54) 

mu 

For u invariants (|53|), ( |5^ coincide with free motion invariants (0), (|3^), 
respectively. 

The propagator of the harmonic oscillator satisfies the system of equations 

dG{q,q',t) . . , .^ ■^dG{q,q',t) 

— tncosuji \- muqsmujt G[q,q,t)=tn , (55) 

oq oq' 

in . dG{q,q',t) . r-i , .^ i^, / ,n (.p^ 

smu;t hg cos cut G[q,q ,t) = q G{q,q ,t), (ob) 

mu oq 

.^ dG{q,q',t) h' ^G{q^ moo' 

dt d'q +^^G{q,q,t), t > 0. (57) 

Integrating Eq. ( ^6]) we find g-dependence of the propagator 



G{q,q',t) = f{q',t)exp{^{q'cotujt - ^)}, (58) 

2h smuit 

in which the factor f{q', t) has to be found using Eqs. ( |55D and (p^). 
Using Eq. ( ^5]) we obtain the following relation 

df{q',t) ,imuj 

—Q^, — = q^— cot ujt f{q,t). (59) 

From this equation g'-dependence of the propagator is found and we have 

G(g,g',t) = iV(t)exp{^[(g2 + g'2)cota;t - ^]}. (60) 

zti sm uit 

Due to Eq. ( |57| ) the time-dependent factor A^( t) satisfies the equation 

N{t) = cot out N{t), (61) 



which gives 

N{t) = ^£=. (62) 



The normahzation constant C may be found from the initial condition ([16| ) using 
Eq. (|50|). So we have the propagator of the harmonic oscillator 

^/ / N I Amuo,, 9 2qq' , , 

V Imn^muot 2n smcut 

Equation (0) for the propagator of the oscillator is the equation in which variables 
q, q', t are not separable. The advantage of the metod of integrals of motion is 
the possibility to use the procedure analogous to the procedure when we solve the 
equation with separable variables. Thus we find the dependence of the propagator on 
different coordinates step by step as in the method of separable coordinates. 



8 Coherent States 

Let us introduce the annihilation and creation operators 

1 , Imuj i . / ,x 

-q+^^p), (64) 



a/2 V h y/hmu 



They satisfy commutation relations 

[a,a]] = i. (66) 
The Hamiltonian of the harmonic oscillator can be rewritten as 

H = huj{a]a+^). (67) 
The integrals of motion linear in the creation and annihilation operators are fl]] , |@] , 

A{t) = e^^*a, (68) 

A]{t) = e-^^*at. (69) 

The operator A{t) is the linear combination of invariants po(t) (p3D and go(^) ( ^41) 
of the form 

= + -T^Poit)]- (70) 

V2 V n yfimuj 

Let us find normalized solutions to the eigenvalue equation 

ai^a{q) = ai^aiq), (71) 



where a is a complex number. Equation ( [TTD rewritten as the differential equation 



yh/muj dq 

has the normalized solution of the Gaussian form 



M,) = (j^Y" exp(-^ - —f + ^—c - -). (73) 

For a = 0, we have the wave function of the ground state 

Mq) =< 9|o >= (7— ) ^^p("^y )• (^4) 

The state \a > which is the normalized eigenstate of the annihilation operator a 

a\a >= a\a > (75) 

is called the coherent state and in coordinate representation the wave function of this 
state is given by Eq. ([73|). Orthogonal number states |n >, n= 0, 1, 2, . . . which 
are eigenstates of the hermitian operator afa, 

afaln >= n\n >, (76) 

satisfy the completeness condition 

00 

\n >< n\ = i. (77) 

n=0 

The coherent state is decomposed as 



lap ^ a" 



|a>=exp ^2^^F>- (78) 



Since 

a\n > = \^\n — 1 >, 



a\\n> = VnTT\n + l>, (79) 



one can check that the series (^) satisfies the eigenvalue equation (ffSl). Due to the 
orthogonality property of the number state basis 

< n\m >= 6nm, (80) 

one has the scalar product of coherent states or matrix elements of unity operator in 
the coherent state basis 

< a\(3 >= exp(-^ ~ ^ + "l^l-^ > ■ 



Using the integral 

J exp(— ax + ox) ax = ^ — exp —, (82) 
and differentiating this equahty as the function of variable a we have 



/oo a/TT 6 

x^ exp(-ax^ + bx) dx = — a^^/^fl + — ) exp — . (83) 
-oo 2 2a 4 



Differentiating equahty (|82D as the function of variable b we have 

/ xexp(— ax + ox) ax = — w — exp— . (84) 

j-oo 2a V a 4a 

The mean value of the oscillator coordinate in the coherent state may be calculated 
using the integral (|84[) 



/oo / 2ri 
Q'lV'a(5')P dq = \ Re a. (85) 
-oo V TTT'^^ 

The mean value of momentum in the coherent state |q; > is also given by the integral 
of the form (|84D , and we have 

/oo Olpaiq^ 
i'ai^) — K — — V2hrnLU Im a. (86) 
-oo oq 

The dispersion of the coordinate in the coherent state 

ag =< a\q'^\a > — (< a\q\a >Y (87) 
may be calculated using the integral (|^) and one obtaines 

h 



2mu! 

The dispersion of the momentum in the coherent state 



ap =< a\p'^\a > -(< a\p\a >)^ (89) 



'p 

may be also obtained by means of the integral (^3]) and one has 

hmuj 

ap = ^— . (90) 

Thus dispersions of momentum and position which in quantum optics are called 
dispersions of quadratures do not depend on the complex number a. The covariance 
of the position and momentum in the coherent state 



'^gp ~ 2 ^ (^\ilP + > — < a\q\a >< a\p\a > (91) 



turns out to be equal to zero 

= 0. (92) 

The decomposition of the coherent state ( ffSD means that the photon distribution 
function of the coherent hght 



I |2ri 

<n|« > P = e-l"l'^ (93) 



is the Poisson distribution function. Thus in coherent states uncertainties of quadra- 
tures are equal to uncertainties in the vacuum state |0 >, and dispersions minimize 
the Heisenberg uncertainty relation 

CTpCT, > — , (94) 

which for coherent state becomes the equality 

(Tp(^q = ^- (95) 

There are no quadrature correlations in the coherent state. The mean value of the 
number of photons in the coherent state is 

oo 

< n >= ^ nP{n) = |a|^ (96) 

n=0 

and the photon number dispersion equals 

cr„ =< > - < n >^= (97) 

since for the Poisson distribution the variance is equal to mean value of a random 
variable. Using the generating function for the Poisson distribution 

oo 

G{z) = ;z"P(n) = exp[{z - l)\a\% (98) 

n=0 

one can easily check these properties due to relations 

dG . , 

(2 = 1) = <n>, 

az 

^-^(z = 1) = <n>'^ - <n> . (99) 
d^z 

The coherent state may be created from the vacuum state by means of the unitary 
displacement operator 

|a >= L)(a)|0 >, (100) 



where 

lap 

D{a) = exp(aat — a*a) = exp( —) expi^aa]) exp(— a*a). (101) 

For matrices A and B of finite and infinite dimensions, there exists the formula 

e^Ae-^ = A+[B,A] + ^[B, [B, A]] + ^[B, [B, [B, A]]] + ... (102) 

In the case 

[A [A,B]] = [B,[A,B]] = 0, 

the following rule holds 

e"" = exp{A + B + ^[A,B]). (103) 

This rule was used in formula ( |101|) since the commutator of the creation and annihi- 
lation operators commutes with these operators. Applying the displacement operator 
to vacuum state and using the equality 

exp(-a*a)|0 >= |0 > (104) 

we reproduce the series (^) since the number state |n > is given by the relation 

in 

|n>=^|0>. (105) 



From the decomposition (^) it is obvious that the vector exp(|ap/2) \a > is an- 
alytical function of the complex variable a, and it is the generating function for 
the number state vector \n >. Using this remark and the definition of Hermite 
polynomials through the generating function 



e 



n=0 



one can obtain the expression for the wave function of the number state ipn{(l) = < 
q\n > comparing the coherent state wave function (^) with the generating function 
(TO- We obtain 

V2"n! V h 

where V'o(q') is given by Eq. (^). 

Using the completeness relation of number states (^) one can obtain the com- 
pleteness relation of coherent states 



\a >< a\ dai da2 = 1, a = ai + ia2, (108) 



since in polar coordinates the integral in the left hand side of equality (|108|) is reduced 
easily to the sum in the left hand side of Eq. ([77|) . 



9 Evolution of Coherent State 



Since the Green function is the kernel of the evolution operator the evolution of the 
coherent state wave function is given by the integral 



M<l,t)= / G{q,q',t)MQ')dq', (109) 
where we use the propagator . The corresponding integral is of the Gaussian form 



(p^) and we have 

Mq,t) = e-^-*/2^~(g), a = ae-'^'. (110) 

Thus in the process of evolution the initial coherent state of the harmonic oscillator 
preserves the property to be the eigenvector of the annihilation operator, but its 
eigenvalue has the time-dependent phase factor corresponding to the classical motion 
of the oscillator in its phase space due to formulae (§^), (§^). It means that in 
the process of evolution dispersions of quadratures do not depend on time and no 
quadrature correlation appears due to the evolution. Having this result it is easy to 
obtain the oscillator propagator in the coherent state basis. We have by definition 

G{a*,(3,t)=<a\U{t)\(3> = 6"^^*/^ < a|/3e-*^* > 

= e-'"*/2exp[-^-^ + a*/3e-'"*]. (Ill) 

Using the property of coherent states to be the generating function for number states 
and expanding the function exp(-'^ — h ^) G{ a*, /3, t) into the power series in 
variables a*, (3 we obtain the obvious expression for the oscillator propagator in the 
Fock basis 

G'(n,m,t) = 5„„e-*-*("+i/^). (112) 



10 Squeezing in Parametric Oscillator 

For the parametric oscillator with the Hamiltonian 

where we take h = m = uj{Q) = 1, there exists the time-dependent integral of 
motion found in [^] 

A = ^[e{t)p - e{t)A. (114) 

where 

e{t) + uj'^{t)e{t) = 0, e{0) = 1, ^(0) = i, (115) 



satisfying the commutation relation 



[A, A]] = 1. (116) 

It is easy to show that packet solutions of the Schrodinger equation may be introduced 
and interpreted as coherent states [Q, since they are eigenstates of the operator A 
(|TT|), of the form 

^„(x,t) = ^^o(x,t)exp{-^^ _^ + ^_}, (117) 

where 

v^o(x,t) = 7r-V4,(t)-i/2exp|^ (118) 

is analog of the ground state of the oscillator and a is a complex number. Variances 
of the position and momentum of the parametric oscillator in the state ( |118|) are 



and the correlation coefficient r of the position and momentum has the value corre- 



sponding to minimization of the Schrodinger uncertainty relation [|13 



a.a, = (120) 

If (Ja; < 1/2 (cTp < 1/2) we have squeezing in photon quadrature components. 

The analogs of orthogonal and complete system of states which are excited states 
of stationary oscillator are obtained by expansion of ( |117] ) into power series in a. We 
have 

and these squeezed and correlated number states are eigenstates of invariant A'^A. 

The function ( |118| ) describes the squeezed vacuum state. The photon distribution 
function in the squeezed vacuum is expressed in terms of the overlap integral 

W{n) = \C{n)\^ = I / ro{x,t)Mx) dx\', (122) 



where the function ipni^) is the wave function of the ra-photon state ( |10?| ). The 
amplitude C{n) can be calculated using the overlap integral 

0(a) = / i:*o{x,t)e^"^ ^^iJa,{x) dx, (123) 



which is generating function for the amphtudes C{n) 



'^(«) = E -F^C^^ (124) 



since according to Eq. (^) the wave function of the coherent state e^"'^^^'^tpa{x) ( [75D 
is the generating function for the n-photon state. The Gaussian integral ( |123| ) is (we 
have h = m = u = 1) 

exp[-- - -(1 + —)x' + v^H, (125) 



and using formula (|82D we obtain 

0^ = (7(0) exp a V, (126) 

where 

m^f^, (127) 

V £ + IS* 2[e* + ze*) 

The probability to have no photons is 

WiQ) = ^ . (128) 

V|2+ |e|2 + 2 



Expanding Eq. ( |125| ) into power series in a we obtain the result that the probability 
to have the odd number of photons n = 2m + 1 is equal to zero 

W{2m+1) = 0, (129) 

and the probability to have even number of photons n = 2m equals 



W{2m) = WiO)p^\fi\^'^. (130) 
[ml Y 



The partial cases of parametric oscillator are free motion ( uj{t) = ), starionary 
harmonic oscillator ( Ci;^(t) = 1 ), and repulsive oscillator ( co'^(t) = —1 ). The 
obtained above solutions are described by the function e{t) which is equal to e{t) = 
1 + it, for free particle, e(t) = e**, for usual oscillator, and e{t) = cosht + isinht, for 
repulsive oscillator. 

If one introduces the squeezing parameter r (for zero correlation) according to 

l^l = e-'\ \e\ = e\ (131) 
the distribution ( |130| ) may be rewritten as 



Another normalized solution to the Schrodinger equation 



^o^mix, t) = 2N^^o{x, t) exp{-^^ cosh — (133) 

where 

iV„. ^ (134, 
2A/cosh laP 



is the even coherent state [|T^ (the Schrodinger cat male state). The odd coherent 
state of the parametric oscillator (Schrodinger cat female state) 

^o.f{x,t) = 2Nf^o{x,t)exp{-'-^ AA_}sinh^^, (135) 

where 

exp(|a|V2) , , 

Nf = ' ' \ (136) 

2A/sinh laP 



satisfies the Schrodinger equation and is the eigenstate of the integral of motion 
(as well as the even coherent state) with the eigenvalue a^. These states are one- 



mode examples of squeezed and correlated Schrodinger cat states constructed in [15 



11 Wigner and Q— functions 



Let us discuss now other possible descriptions of the quantum state. There exist 
two important functions which we describe for multimode system with N degrees of 
freedom. The Wigner function of a system Vr(p, q) = Vr(Q) is expressed in terms 
of density matrix in dimensionless coordinate representation as (see, |1l6| ) 



W{p, q) = y p(q + ^, q - ^) exp(-zpu) du 
The inverse transform is 

1 



p(x, x') 



(2vr: 



X + x' 

W{ — - — , p) exp[ip(x — x')] dp. 



(137) 



(138) 



The Q-function is expressed in terms of Wigner function through the 3N- 
dimensional integral transform 



Q(B) 



i2n 



X + x' 

$b(x, x', p)W{ , p) (ix (ix' dp, 



(139) 



with the kernel 



$b(x, x', p) = tt exp[ip(x — x' 



1 

-B((Tjv. + l2Jv)B - ^ + v^Ba^.X], (140) 



where the 2N-vector X = (x, x') is introduced. If one has the Q-function the Wigner 
function is given by the integral transform 



I r ^ ~ ~ ~ 

W^(p,q) = ^ iUd^l^k Afe duk^k{uk,B)}Q{B), (141) 
k=i 

where the argument of the Q-function B is replaced by the 2N-vector with complex 
components 

B = ( A, /?2, • • • , f^N, 7l, 72, • • • ) 77v), 
and the kernel has the form 

Muk,B) = 7r-'/\xp[-\(3,\'-\^,\' + V2{q, + ^)(3, + V2{q,-^)jl 



- ^fe + Y)^-^fe-Y)^-f -|^-W + 7,/5-]. (142) 



The density matrix in coordinate representation is related to the Q-function 

p{^,^')=n-''' J{Y[d'(3, rfS0fc(B)exp[--(4 + xL^)]}Q(B), (143) 

k=l 

where the kernel of the transform is 

0,(B) = 7r-V2 exp[-|/?,|2 - |7,|2 + V2xA + ^^x'^ll " f - ^ + 7^/5.1- (144) 

The evolution of the Wigner function and Q-function for systems with quadratic 
Hamiltonians for any state is given by the following prescription. Given the Wigner 
function iy(p, q, t = 0) for the initial moment of time t = 0. Then the Wigner 
function for time t is obtained by the replacement 

W^(p,q,t)= l^(p(t), q(t), t = 0), 

where time-dependent arguments are linear integrals of motion of the quadratic sys- 
tem found in [m, and |jl8[. The same ansatz is used for the Q-function. Namely, 
given the Q-function of the quadratic system Q(B, (t = 0)) for the initial moment 
of time t = 0. Then the Q-function for time t is given by the replacement 

Q(B, t)= g(B(t), t = 0), 

where the 2N- vector B(t) is the integral of motion linear in the annihilation and 
creation operators found in and JTSf. This ansatz follows from the statement 
that the density operator of the Hamiltonian system is the integral of motion, and 
its matrix elements in any basis must depend on appropriate integrals of motion. 
In particular, the Wigner function and Q-function depend just on linear invariants 



found in Hi, HI and lT8f 



As an example we consider the Wigner function of the oscillator ground state with 
the density matrix in coordinate representation (we reconstruct below the dimension- 
ality of variables) 

Po{x,x') = ^— exp[-- — (x^ + x'^)]. (145) 
Using the definition of the Wigner function we have 

*^»'''' = V r7 - 2' + + 2*1 - X> 

and calculating the Gaussian integral we obtain 

W.M^2.M-£^-'l^). (147) 



The Husimi Q-function equals 



= 1 +,^). (148) 

Since Q-function is the diagonal matrix element of the density operator in coherent 
state basis it is easy to obtain the Q-function of the thermal equilibrium state of the 
oscillator at temperature T with density operator 

p(T) = 2sinh^exp[-^(ata + 1)]. (149) 

The operator Zp{T), where Z is the partition function 

Z = (2sinh^)-\ (150) 

coincides with the evolution operator U{t) if one takes it = h/T. This equality 
gives in coherent state basis the Q-function of the oscillator at temperature T 

Q.(,.p) = 2smh-exp{-- - -(£- + - exp(--)l}, (151) 

if one uses formula ( p.ll| ) for the oscillator propagator in coherent state basis 



\3 >= IqJ—- + I , > . 



12 Multimode Mixed Correlated Light 

The most general mixed squeezed state of the N-mode hght with a Gaussian density 
operator p is described by the Wigner function W{p, q) of the generic Gaussian 
form. 



iy(p,q) = (detM)"2 exp 



1 



(Q- < Q >)M i(Q- < Q >) 



(152) 



where 2N-dimensional vector Q = (p, q) consists of N components pi, . . . , and 
N components qi, Qn] operators p and q being quadrature components of 

the photon creation af and annihilation a operators (we use dimensionless variables 
and assume ^ = 1): 



P 

q 



at 



iV2 ' 
a + af 



V2 



2N parameters < Pi > and < qi >, i = 1,2, 
are average values of quadratures, 



(153) 



.N, combined into vector < Q>, 



< P > 

< q > 



Tr pp, 
Tr pq. 



A real symmetric dispersion matrix M consists of 2N^+N variances 
1 



M 



QaQp + QpQa) — \Qa) \Q 0) , 



a,P = l,2, 



2iV. 



(154) 



(155) 



They obey certain constraints, which are nothing but generalized uncertainty relations 
i- 

The photon distribution function of this state 



Vn = Tr p|n >< n| 



n = (?2i,n2, . . .,12^), 



(156) 



where the state |n > is photon number state, which was calculated in ||T9|, |22| and 
it is 

= Vo^^^^. (157) 



n! 



The function Hj^^{y) is multidimensional Hermite polynomial. The probability to 
have no photons is 



det ( M + -I2N 



exp 



- < Q > (2M + I 



2N 



< Q > 



(158) 



where we introduced the matrix 

R = 2U^(l + 2M)-iU*-a7Vx, (159) 

and the matrix 

In 



^iv.= ° ^n- (160) 



The argument of Hermite poljTiomial is 

y = 2U*(l2^ - 2M)-i < Q >, (161) 
and the 2N-dimensional unitary matrix 

U=^f-;^^ (162) 

is introduced, in which I^r is the NxN identity matrix. Also we use the notation 

n! = ni\n2\ ■ ■ - riN^.. 

The mean photon number for j-th mode is expressed in terms of photon quadrature 
means and dispersions 

< >= ^(^P,P, + ^9,% -^) + l{< Pj >' + < >)• (163) 

We introduce a complex 2N-vector B = (/3i, /?2, • • • , Pn, P*, ■ ■ ■ , Pn)- Then 
the Q-function Jl^ is the diagonal matrix element of the density operator in coherent 
state basis \ Pi, P2, ■ ■ ■ , Pn > ■ This function is the generating function for matrix 
elements of the density operator in the Fock basis |n> which has been calculated in 
In notations corresponding to the Wigner function ( 152 ) the Q-function is 



Q(B) = Vo exp[-^B(i? + a^,)B + Bi?y]. (164) 



Thus, if the Wigner function (|152D is given one has the Q-function. Also, if one has 



the Q-function ( |164| ), i.e. the matrix R and the vector y, the Wigner function may 
be obtained due to relations 

M = U*(i? + a^,)"'Ut-i, 
<Q> = U*[l-(i? + a^,)-V;vJy. (165) 
For pure squeezed and correlated state with the wave function 

\Ef = A^exp[— xmx + cx], (166) 



where 

N= [det(m + m*)]i/V-^/^exp{^(c + c*)(m + m*)-i(c + c*)}, (167) 

the symmetric 2Nx2N-matrix R determining Q-function has the block-diagonal 
form 

^=(0 r*)' ^^^^^ 
The NxN-matrix r is expressed in terms of symmetric matrix m 

r* = 1 - (m + l/2)~\ (169) 

and the 2N-vector y = (Y, Y*) is given by the relation 

Y* = i=(m-l/2)-ic. (170) 

Corresponding blocks of the dispersion matrix 

M=( "^P^ \ (171) 



are 



o-qq = ^(m + m*)"\ 

Cpq = -(m — m*)(m + m*)~"^. (172) 



The probability to have no photons is 

^ _ [det(m + m*)]^/2 
° ~ |det(m + l/2)| 

® exp{^(c + c*) (m + m*)-^{c + c*) + ^[c(m + l/2)"^c + c*(m* + l/2)"^c*]}. (173) 

Multivariable Hermite polynomials describe the photon distribution function for the 
multimode mixed and pure correlated light |T^, ||2T|, p2|. The nonclassical state of 



light may be created due to nonstationary Casimir effect and the multimode 
oscillator is the model to describe the behaviour of squeezed and correlated photons. 



13 Multivariable Hermite Polynomials 



For parametric forced oscillator, the transition amplitude between its energy levels 
has been calculated as overlap integral of two generic Hermite polynomials with a 
Gaussian function (Frank-Condon factor) and expressed in terms of Hermite polyno- 
mials of two variables [§]. For N-mode parametric oscillator, the analogous amplitude 
has been expressed in terms of Hermite polynomials of 2N variables, i.e. the overlap 
integral of two generic Hermite polynomials of N variables with a Gaussian function 



(Frank-Condon factor for a polyatomic molecule) has been evaluated in 0, [jl8|. The 
corresponding result uses the formula 

[n = ni,n2, . . . ,nN, m = mi,m2, . . .ruN, m^, ra, = 0, 1, . . .) 
r 7r^/2 1 

/ H^^\^)H^\A^+d) exp(-xmx+cx)rfx = exp{-cm-^c)H^l{y), (174) 

J V det m 4 

where the symmetric 2Nx2N-matrix 

with NxN-blocks R2, R12 is expressed in terms of symmetric NxN-matrices 
i?, r, m and NxN-matrix A in the form 

Ri = R- -Rm-^R, 
2 

R2 = r — —rAm^^Ar, 

R12 = -Katu^^R. (176) 

Here the matrix A is transposed matrix A, and Ru is transposed matrix R12. The 
2N-vector y is expressed in terms of N-vectors c and d in the form 



72 



where N-vectors yi and y2 are 



yi = ^{Rfn ^ + m ^R)c 
1 

y2 = -{rAm~^ + m'^^Ar)c + rd. (178) 



For matrices R = 2, r = 2, the above formula (|174D yields 

N N 



J {f[ HnXxi)HjnXYl ^ikXk + c?i)} exp(-xmx + cx)dx 



i=l k=l 



77^/2 1 

= ^==exp(-cm-ic)i7i^>(y), (179) 

with NxN-blocks Ri, R2, Ru expressed in terms of NxN-matrices m and A in 
the form 

Ri = 2(1 -m-^), 
R2 = 2(1-Am-^A), 

R12 = -2Am-\ (180) 
The 2N-vector y is expressed in terms of N-vectors c and d in the form ( |1 77| ) with 

yi = m'^c, 

ys = ^(Am-^ + m~^A)c + 2d. (181) 
If the symmetric matrix p has the block-diagonal structure 



P 



Ri 
R2 



with the symmetric SxS-matrix Ri and the symmetric (2N-S)x (2N-S)-matrix R2, 
the multivariable Hermite polynomial is represented as the product of two Hermite 
polynomials depending on S and 2N-S variables, respectively, 

HfHy) = HifHy,)Hif^l{y2), 

where the 2N-vector y has vector-components 

y = ( yi, y2), 

and the 2N-vector k has components 

k = (ns, n2Ar_s) = ( "-i, • • • , "-5, ns+i, n2N)- 

The partial case of this relation is the relation for Hermite polynomials with the 
matrix R with complex conjugate blocks Ri = r, R2 = r*, and complex conjugate 
vector-components yi = yg 

HtHy) = \HifHyi)\', s = N. 

The calculated integrals are important to evaluate the Green function or density ma- 
trix for systems with quadratic Hamiltonians. Partial cases of multivariable Hermite 
polynomials determine some other special functions [BTI], P2 . 



14 Multimode Even and Odd Coherent States 

We define multimode even and odd coherent states (Schrodinger cat male states and 
Schrodinger cat female states, respectively) as 



I A± >= A^±(| A > ± I -A >), (182) 

where the multimode coherent state | A > is 

I A >=| «!, 0:2, . . . , a„ >= -D(A) I >, (183) 

and the multimode coherent state is created from the multimode vacuum state | > 
by the multimode displacement operator -D(A). The definition of multimode even 
and odd coherent states is the obvious generalization of the one-mode even and odd 
coherent state given in 0, [Q. Normalization constants for multimode even and odd 
coherent states become 

plA|V2 



2a /cosh I A p 



elAp/2 

iV_ = , ^ (184) 
2A/sinh A 



|2 

where A = ( ai, 0^2, . . . , an) is a complex vector and its modulus is 

n 

\A\^=J2\am\^. (185) 

m=l 

Such multimode even and odd coherent states can be decomposed into multimode 
number states as 

I A± >= N± "1 ^ ^_^^„,+n2+...n„) I ^ >^ ^^gg^ 

n V ^1 • ■ ■ ■ V TT-n- 

where | n >= \ rii, n2, . . . , n„ > is multimode number state. Also from Eq. ( |182| ), 
we can derive an important relation for the multimode even and odd coherent states, 
namely. 



tti I A_i_ > = aa/tanh | A ^ | A_ >, 



tti I A_ > = ajycoth I A P I A+ > . (187) 

The probability of finding n photons in multimode even and odd coherent states can 
be worked out with the help of Eq. (TE^) 



-^+(n) = 7 — ^T7 — 7 — , , , ,0 , ni + n2 H h n„ = 2/c, 



ai 


2ni 


a2 


2*^2 .. . 


an 


2n„ 


(ni!)(n2! 
1 ai 


)■■■ 

^2 


(nj.) cosh 

2«2 . . . 1 ^ 

1 


A |2 

2n„ 



P-(n) = , ' ' \ ,^ .' I ' I, , ni + n2 + --- + n„ = 2A: + l. (188) 

(ni!)(?22') ■ ■ ■ (n„!) smh | A p 



Multimode coherent states are the product of independent coherent states of each 
mode, and photon distribution function is the product of independent Poissonian 
distribution functions. But in the present case of muhimode even and odd coherent 
states we cannot factorize their multimode photon distribution functions due to the 
presence of the nonfactorizable cosh | A p and sinh | A p. This fact imphes the 
phenomenon of statistical dependences of different modes of these states on each 
other. 

In order to describe properties of distribution functions from Eq. ( |188| ) we will 
calculate the symmetric 2Nx2N dispersion matrix for multimode field quadrature 
components. For even and odd coherent states, we have 

< A± I ttittk I A± >= Oittfe, (189) 

and complex conjugate values of the above equation for < A-t | a\a\ \ A± >. Since 
the quantity < A± | | A± > is equal to zero the above equation gives two NxN 
blocks of the dispersion matrix. For other two NxN blocks of this matrix, we have 
for multimode even coherent states 

ol'i =< A+ I \{a\ak + Ofcaf) | A+ >= a*afctanh | A |^ +\Sik, (190) 

and for multimode odd coherent states 

(T"t =< A_ I \{a\ak + akol) \ A >= a*ak coth | A |^ +^Sik. (191) 

(a\ak) '2 2 

For the dispersion matrix, the mean values of the photon numbers rii = a|aj for 
multimode even and odd coherent states are the following 



< A+ 


\ni\ A+> = 


p tanh 


1 Ap, 




< A_ 


\ Uil A_ > = 


a, P coth 


AT- 


(192) 



Taking into account the above equation the symmetric NxN dispersion matrices for 
photon number operators can be obtained from the above given distribution functions 
for multimode even and odd coherent states. By defining 

=< A± I UiUk I A± >, (193) 

corresponding expressions in such states are 

^ifc — I '^i n '^k P sech ^ I A p + I p tanh | A ^ 6ik, 

= - I ai pl ak P csch ^ | A ^ + | ai ^ coth | A ^ 6ik. (194) 

As nondiagonal matrix elements of the dispersion density matrix are not equal to zero, 
we can predict that different modes of these states are correlated with each other. In 



other words, as we have mentioned before, there exist some statistical dependences 
of different modes on each other. 

Another interesting property for muhimode even and odd coherent states is the 
Q-function, and it can be obtained in the following manner. First of all the density 
matrices for multimode even and odd coherent states are 

P±=|A±><A±|, (195) 

then the Q-function can be calculated as 

Q+(B,B*) = <B|p+|B> 

= 4Ar2e-(l^l'+l^l') I cosh(AB*) |^ 

Q_(B,B*) = < B I p_ I B > 

= 4Ar!e-(l^l'+l^l') I sinh(AB*) \\ (196) 

where | B >=| Pi, P2, • • • , /?n > is another multimode coherent state with multimode 
eigenvalue B = (/5i, /?2, • • • , Pn)- We call these functions for even and odd coherent 
states as the Q-functions for Schrodinger cat states. The Q-function for single-mode 
odd coherent state shows the crater type behaviour for small values of the quantity 
I a I and for its larger values the Q-function begins to split into two peaks in a similar 
manner as in case of even coherent states P3 . 



The Wigner function for multimode coherent states is |]l| 

I A P I B P 

Wa,b = 2^ exp[-2ZZ* + 2AZ* + 2B*Z - AB* - ^ ^ - ^ (197) 



where 

Z = ^. 

For even and odd coherent states, the Wigner function is 



(198) 



W^A±(q, p) = I A^i P [Vr(A,B=A)(q, p) ± Vr(A,B=-A)(q, p) 

± W^(-A,B=A)(q,p) + W^(-A,B=-A)(q,p)], (199) 



where explicit forms of N± are given in Eq. ( |184| ). For multimode case, we use 
following notations 



AZ* = aiZl + 0:2^2 + ■ ■ ■ cunZ*, 

ZZ* = ZiZ; + Z2Z; + ■ ■ ■ + ZnZ:,. (200) 

The photon distribution function gives the probability of finding 2k photons for 
two-mode even coherent state, and is defined as 



(1 "1 


2 + 1 


^2 


2yk 


(2A;)!cosh(| ai 


P + 


Ci2 





P+W = .O.N 'An '1/ . i.N ^ (201) 



where 2k = ni + n2, for both rii and n2 to be even or odd numbers. Similarly for 
two-mode odd coherent state it gives the probability of finding 2k + 1 photons 



(1 "1 


2 + 


a2 








(2A; + l)!sinh(| ai 


2 + 


a2 





P-(2fc + l)= / . , I. (202) 



For this case, rii is even and n2 is odd number, or vice versa. For single-mode case, 
the photon distribution functions demonstrate super and sub-Poissonian properties 
for even and odd coherent states, respectively. The same conclusion may be drawn 
for two-mode (and multimode) even and odd coherent states. 

15 Generic Quadratic Systems 

The generic mechanical or optical linear system with energy operator has the Hamil- 
tonian 

H = -QB{t)Q + C{t)Q, (203) 

where we use 2N-vectors Q = (pi, P2, ■ ■ ■ ,Pn, Qi, 12, ■ ■ ■ , Qn), and C(t) as well 
as, 2Nx2N-matrix B{t), the Planck constant h = 1. This system has 2N linear 
integrals of motion which are written in vector form 

Qo(t) = A(t)Q + A(t). (204) 

The real symplectic matrix A{t) is the solution of the system of equations 

A{t) = A{t)^B{t), 

A(0) = 1, (205) 

where the real antisymmetrical matrix S is 2N-dimensional analog of the Pauli 
matrix iay, and the vector A(t) is solution of the system of equations 

A(t) = A(i)EC(t), 

A(0) = 0. (206) 
If for time t = 0, one has the initial Wigner function of the system in the form 

W{p,q,t = 0) = Wo{Q), (207) 
the Wigner function of the system at time t is 

W{p,q,t) = Wom)Q + A(t)]. (208) 



The Hamiltonian (g03D may be rewritten in terms of creation and annihilation oper- 
ators 

H = -AD{t)A + E{t)A, (209) 



where we use 2N-vectors A = ( ai, 02, . . . , a^, ctif, ci2t, • • • , «Aft), and E(t) as 
well as, 2Nx2N-matrix D{t). This system has 2N-linear integrals of motion p[ 
which are written in vector form 

Ao(t) = M(t)A + N(t). (210) 

The complex matrix M[t) is the solution of the system of equations 

M{t) = M{t)aD{t), 

M(0) = 1, (211) 

where the imaginary antisymmetric matrix a is 2Nx2N-analog of the Pauli matrix 
— (Ty, and the vector N(t) is solution of the system of equations 

N(t) = M{t)aE{t), 

N(0) = 0. (212) 
If for time t = 0, one has the initial Q-function of the system in the form 

Q{a,a*,t = 0) = Qo{a), (213) 
the Q-function of the system at time t is 

Q{a,a*,t) = Qo[M{t)a + -N{t)]. (214) 

Here a = ( q + ip) / 

For time-independent Hamiltonian (|203|) , the matrix A{t) is 

A(t) = exp(Sfit), (215) 

and the vector A(t) is 

A(t) = f'expiEBr) E C(r) dr. (216) 
Jo 

For time-independent Hamiltonian ( |209|) , the matrix M{t) is 

M{t) = exp{aDt), (217) 

and the vector N{t) is 

N(t) = ^ expicxDr) a E(r) dr. (218) 
Jo 

For time-dependent linear systems, the Wigner function of generic squeezed and 
correlated state has Gaussian form and it was calculated in iTH. 



16 Optical Homodyne Tomography 



In it was shown a relation of marginal distribution w{X, B) for homodyne output 
variable 

X(e) = gcosG -psinG, (219) 
and the Wigner function W{q,p) of the form 



w{X, Q) = J W{x cos 9 — f sin 9, X sin 9 + f cos 9) dv, 



where 



/oo 
w{x,e) dx = 1. 
-oo 

The inverse relation which is Radon tomographic transform, 
Wiq,p,s^O) -- ^ 



(220) 



(221) 



, dx dy dQ w(X, 9) 

4:71^ J-OO Jo 



X exp[(s|/78) + z?/(X - gcos9 -psin9)] |?/|, (222) 



gives the possibility measuring the homodyne variables (|219|) to measure the Wigner 
function of the system. In EE\ the method was extended to measure the Wigner 
function by measuring the marginal distribution w{X, n, u, 6) for the variable 



X = + up + S, 
since it is easy to show that the Wigner function is 

W{q,p) = (27r)^s^ exp(isX) w{X, sq, sp, s). 



(223) 



(224) 



The function w{X, 171,11, s) is the Fourier transform of the marginal distribution of 
w{X,fi, u, 6). 



17 Conclusion 

The discussed method of the time-dependent integrals of motion may be applied 
to systems with varying parameters, for example, to the model of the creation of 
photons and squeezing phenomenon in frame of nonstationary Casimir effect using the 
generalized harmonic oscillator (see, [27]). The analogous phenomenon of reduction 
of quantum noise for oscillating mirrors has been discussed in The behaviour 



of the system containing moving mirrors and light beams has been discussed in p9 



and using even and odd coherent state light (Schrodinger cat light) in interferometric 
gravitational wave detectors has been suggested in []30| . 

The coherent states [^, which describe classical states of light have the 



Poisson photon statistics. Nonclassical states like squeezed states |0, [^, ||36[ or 



generalized correlated states have nonpoissonian light statistics. Different kinds 
of nonclassical states are studied in 



38 



There is possibility to introduce different types of nonclassical states using different 
types of superposition of complete set either of photon number states or of coherent 
states. All such states which may be constructed experimentally are interesting, 
and many of these states may be created in situations with nonstationarities or in 
situations with nonlinearities. 
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